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Abstract In recent work, we analyzed the dynamics of spherical and non- 
spherical acoustic perturbations of the Michel flow, describing the steady ra¬ 
dial accretion of a relativistic perfect fluid into a nonrotating black hole. We 
showed that such perturbations undergo quasi-normal oscillations and com¬ 
puted the corresponding complex frequencies as a function of the black hole 
mass M and the radius Tc of the sonic horizon. It was found that when Vc is 
much larger than the Schwarzschild radius th = 2GM/c^ of the black hole, 
these frequencies scale like the surface gravity of the analogue black hole as¬ 
sociated with the acoustic metric. 

In this work, we analyze the Newtonian limit of the Michel solution and 
its acoustic perturbations. In this limit, the flow outside the sonic horizon 
reduces to the transonic Bondi flow, and the acoustic metric reduces to the one 
introduced by Unruh in the context of experimental black hole evaporation. 
We show that for the transonic Bondi flow, Unruh’s acoustic metric describes 
an analogue black hole and compute the associated quasi-normal frequencies. 
We prove that they do indeed scale like the surface gravity of the acoustic black 
hole, thus providing an explanation for our previous results in the relativistic 
setting. 
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1 Introduction 

Studying the accretion of matter into a black hole is a highly relevant sub¬ 
ject in general relativity and astrophysics. From a purely theoretical point of 
view, it allows us to understand the dynamical evolution of black holes, a 
topic that has played a prominent role in the development of black hole ther¬ 
modynamics [UiniElIl] and is also directly related to the black hole stability 
problem On the astrophysical front, the study of accretion processes are 
utterly important since they provide a mechanism by which black holes can 
be “observed”. Traditionally, these observations are based on measuring the 
electromagnetic spectrum emitted by the accretion disk and comparing this 
spectrum to theoretical models (see Ref. [7] and references therein). Recently, 
millimeter-wave very-long baseline interferometric arrays such as the Event 
Horizon Telescope [S] have been able to resolve the region around Sagittarius 
A*, the supermassive black hole lying at the center of our galaxy, to scales 
smaller than its gravitational radius [9]. It is expected that the comparison 
between the observations to calculated images of the black hole shadow and 
the sharp photon ring surrounding it lead to tests for the validity of general 
relativity in its strong (albeit stationary) regime. 

Another astrophysical scenario in which accretion comes into play is in a 
black hole-neutron star binary. For supermassive black holes, the star remains 
intact as it spirals around the black hole, and eventually it falls into it as 
a consequence of energy loss through gravitational radiation. However, for 
less massive black holes, more precisely, for black holes whose tidal radius is 
comparable or larger than the radius of its inner stable circular orbit, the star 
is tidally disrupted before it falls into the black hole, and an accretion disk 
forms. The formation of this disk has a significant impact on the gravitational 
wave signal, and so the study of this process is relevant for gravitational wave 
experiments such as LIGO [TD], VIRGO [TT] and KAGRA |T1], see [T3] and 
references therein for a recent review on such binary systems. 

In this article, we focus on the case of nearly spherical and steady-state 
accretion by a nonrotating black hole. Although this simple model does not 
say much about the scenarios described so far, it is still relevant in a number 
of astrophysical problems, like the accretion of interstellar matter by a stellar 
black hole [MllIlllT], or the accretion of dark matter by a supermassive or 
stellar size black hole [Hllil. For recent work about radial accretion in a 
cosmological background, see also [Biiiiiiniiii]- 

In recent work mm, we provided a thorough analysis for the Michel 
flow m, describing the spherical, steady-state accretion into a nonrotating 
black hole. Unlike previous results which mainly focused on the case of a 
Schwarzschild background and a polytropic perfect fluid, our result in [23] 
shows that given an arbitrary static, spherically symmetric black hole back¬ 
ground satisfying certain reasonable conditions, and given a perfect fluid on 
this background with an arbitrary equation of state satisfying certain general 
conditions, there exists for each value rioo of the particle density at infinity 
a unique accretion flow which is everywhere regular at and outside the event 
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horizon. As it turns out, this flow is necessarily transonic, the fluid’s velocity 
(as measured by static observers) being subsonic in the asymptotic region and 
supersonic close to the black hole. 

In an extension of this work, in [24] we studied small, spherical and non- 
spherical acoustic perturbations of the Michel flow and showed that they ex¬ 
hibit quasi-normal acoustic oscillations. Their appearance is most easily un¬ 
derstood by following the method by Moncrief [25] who showed that the dy¬ 
namics of the acoustic perturbations are described in a very elegant manner 
by a wave equation on an effective curved background geometry determined 
by the acoustic (or sound) metric. As it turns out, the geometry described by 
the acoustic metric has precisely the same features as a static and spherically 
symmetric black hole, where for the particular case of Michel accretion the 
horizon is located at the sonic sphere, and thus it describes a natural astro- 
physical analogue black hole. In view of these observations, the appearance of 
quasi-normal acoustic oscillations is no surprise. 

On the other hand, as discussed in |21|, the calculations of the correspond¬ 
ing modes and frequencies is technically challenging, since the acoustic metric 
is not known in explicit closed form. Consequently, traditional methods for 
computing the quasi-normal frequencies, such as Leaver’s method [26] cannot 
easily be applied to this problem, and we had to design a new computation 
method to calculate the modes. Based on this new method, which is described 
in [24] and in more detail in EH, we computed the quasi-normal complex fre¬ 
quencies s = a -\- ioj as a, function of the black hole mass M and the radius 
Tc of the sonic sphere for a polytropic fluid with adiabatic index 7 = 4/3. We 
noticed that for = 2GM/c^, both the decay rate a and the frequency 

w of the oscillatory part of the mode scale like the surface gravity k of the 
acoustic black hole. 

The main goal of this article is to provide an explanation for this behav¬ 
ior by reconsidering the accretion problem in the Newtonian limit. Here, by 
“Newtonian limit” we mean the realistic case in which the speed of sound 
of the fluid Vao at infinity is much smaller than the speed of light c. In this 
limit, as further discussed below, it can be shown that the accretion flow can 
be described to very good accuracy by the Newtonian equations in the region 
outside (and including) the sonic horizon. The radial, steady state Newtonian 
solutions were studied in detail by Bondi |I4| a long time ago, and as it turns 
out the Michel flow reduces to the transonic Bondi flow in the Newtonian 
limit |7]. 

After presenting a very brief, though self-contained derivation of the tran¬ 
sonic Bondi flow solution, we derive the equations describing small acous¬ 
tic perturbations thereof. Interestingly, the nonrelativistic description of this 
problem is still described by a wave equation with an effective curved metric 
background in four spacetime dimensions. The acoustic metric, which can be 
obtained from the Newtonian limit of Moncrief’s sound metric, has been ob¬ 
tained by Unruh |28| in the context of an experimental black hole evaporation 
proposal. We discuss the geometry of the acoustic metric in the Newtonian case 
and compute the quasi-normal frequencies of the associated wave equation. In 
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particular, we prove that these frequencies scale like the surface gravity of 
the acoustic hole and thereby provide an explanation for the results obtained 
in [23] in the relativistic case. Furthermore, we also discuss the eikonal ap¬ 
proximation which connects the quasi-normal frequencies to the properties of 
the unstable circular null geodesics and find that it yields a fairly good de¬ 
scription for the decay rates and the change in the oscillation frequencies as 
the angular momentum increases by one. Finally, we analyze the behavior of 
the quasi-normal frequencies s with respect to the adiabatic index 7 . 


2 Transonic Bondi flows 


The equations of motion describing a nonrelativistic, barotropic fluid in an 
external gravitational potential (j>, consist of the continuity equation 

p + V-(pv) = 0 (1) 


and the Euler equation 


V -1- (v • V)v -I- V [h{p) -I- (^] = 0, 


( 2 ) 


where p, v and h{p) = f dp/p refer to the density, the velocity and the en¬ 
thalpy per unit mass of the fluid, respectively, and where a dot denotes partial 
differentiation with respect to time. For the following, we assume a polytropic 
equation of state, p = Kp'^^ with 7 > 1 the adiabatic index, in which case 
h{p) = Ug /(7 — 1 ), with Vs = \Jdp/dp = the local speed of 

sound in the fluid. 

For the purpose of the present article it is sufficient to consider the case 
where the fluid is irrotational, such that v = Vip with ip the fluid potential. 
In this case, Eq. (§ can be integrated 

V'+ +/i(p) + (/i = const. (3) 


This equation may be used to express p in terms of first derivates of ip, and 
when introduced into Eq. 0 this yields a nonlinear wave equation for ip. For 
a mathematical discussion regarding the formation of shocks based on the 
relativistic generalization of this wave equation, see [Mlisn]. 

For the particular case of a radial fluid flow accreted by a central body of 
mass M, Eqs. 01 ) can be integrated and give 


Airr^pu = const. 


1 2 
XM + 


GM 


7-1 


= const. 


( 4 ) 

( 5 ) 


where u = v'" is the radial component of the velocity and G is Newton’s 
constant. The constant in the first equation is the mass accretion rate M. 
Eqs. (4[5 1 describe an algebraic system of two equations for the three unknowns 
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{r,p,u). In order to analyze this system, following Refs. we introduce 

the dimensionless quantities 


u 




2 .^ 7-1 ' M 

■ 7-hi’ ■ 47r(GM)2poo’ 


where poo and Voo refer to the density and sound speed, respectively, at infinity. 
With this notation Eqs. Q can be rewritten as 

= A, ( 6 ) 


and Eq. 0 yields 


f{„) = A 


(7) 


with the functions / and g defined by 


m := . 

S(x) := x*' (i + . 

The function / : (0,oo) —?> K diverges to -too as 0 or j/ —oo and has 
a global minimum &t u = = 1. Provided that 1 < 7 < 5/3, the function 
g : ( 0 ,oo) —)• K also diverges to -hoo as a: —)■ 0 or a; —>■ 00 , and in this case g 
has a global minimum at a; = Xc = (5 — 87 )/4. 

Although there are many other accretion flow solutions, in this article we 
are only interested in transonic accretion flows, the flow being subsonic {v < 1 ) 
for large x and supersonic {12 > 1) close to the attractor. At the sonic sphere 
1 / = r'c = 1 ; where the flow changes from sub- to supersonic, the derivative 
of / vanishes, and according to Eq. Q the derivative of g also has to vanish 
there for the flow to be regular. Consequently, the transonic flow solution has 
to pass through the point (x, v) = (xc, r'c), and Eq. Q yields 


A = Ae 


1 / 5-37 
4 2 


5-3-7 

2(7-1) 


( 8 ) 


The compression rate at x = Xc is 


— 




(9) 


Outside the sonic sphere, the flow is uniquely determined by Eq. Q where x 
and v are restricted to the intervals (xc, 00 ) and ( 0 , i^c), respectively, on which 
the functions g and / are monotonously increasing. Multiplying both sides of 
Eq. Q by and taking the limit x —>■ 00 , we find that along the flow 


X V —'f A, 
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in this limit, and consequently it follows from Eq. ([^ that z 1 and Vs —> Voc, 
as required. Therefore, the radial velocity of the flow, u = Vsi', converges to 
zero implying that the flow is at rest at infinity. 

When Voo is much smaller than the speed of light c, the transonic Bondi 
flow describes with very good accuracy the accretion flow into a Schwarzschild 
black hole in the region outside the sonic sphere. Indeed, when Vao/c 1, 
and 1 < 7 < 5/3 the sonic sphere is located at a radius Tc satisfying Tc/th = 
(5— 37)/(8z;^/c^) (see for instance [23]), with th = 2GM/c^ the Schwarzschild 
radius of the hole, and hence the exterior of the sonic sphere lies entirely in the 
weak held region. Furthermore, since u/c = {voo/, the magnitude 
of the fluid’s three-velocity, as measured by static observer, is much smaller 
than the speed of light in the region exterior to the sonic sphere. 

The relativistic generalizations of Eqs. are (see, for instance mm) 


hrel (^ 


rn 


= const. 

(10) 

= const. 

(11) 


with hrei {n) the relativistic specihc enthalpy and n the particle density. In the 
nonrelativistic limit, u'' ^ u, n ^ p/m with m the particle mass, hreiin) —^ 
m[c^ -I- h{p)], and using h{p) , rn/r 1, u'^/c <C 1 for all r > Vc one sees 
that Eqs. ( I0|1I ) indeed reduce to Eqs. (4|5|. 


3 Acoustic metric and the geometry it describes 

In this section, we consider spherical and nonspherical linearized acoustic per¬ 
turbations of the Bondi how, keeping the gravitational potential (p = —GM/r 
hxed. The dynamics of these perturbations can be obtained by replacing 

p ^ p + eSp, Ip Ip + eSip 

in Eqs. ( |1|3| ) and linearizing with respect to the small parameter e, where 
(p, V = S/ip) describes the (spherical and steady state) Bondi how solution 
and {6p,6v = Vdip) the (non-spherical and time-dependent) acoustic pertur¬ 
bations. The resulting equations are 

j + ■VS^ + AS^ = 0, ( 12 ) 

Sip + pv ■ \')S'p + Vg—= 0, (13) 

where we have used the identity pdh/dp = v/ and the background equation 
V • (pv) = 0. Eliminating 5p/p from these two equations one obtains the 
following wave equation for Sp) |55|: 


'G\@5'p = 0 , 


( 14 ) 
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where the box operator is computed with respect to the acoustic metric 


0 = — \—vldt^ + 5ij[dx^ — dt){dx^ — v^dt)\ . 


(15) 


The fact that the dynamics of acoustic perturbations of an irrotational New¬ 
tonian fluid flow can be described by a wave equation on an effective curved 
background has been pointed out by Unruh [23 in the context of black hole 
analogue models and the possibility of measuring the analogue of Hawking ra¬ 
diation in these models. For recent experiments realizing such analogue models 
through fluid flows in the laboratory, see for example Ref. m- 

Interestingly though, it seems that the description of the propagation of 
acoustic perturbations through an effective curved geometry had been found 
earlier by Moncrief [25j in an astrophysical context. In his work, Moncrief 
showed the propagation of acoustic linear perturbations of a relativistic, po¬ 
tential flow can be described by a wave equation on a curved background 
described by the sound metric. In the non-relativistic limit, Moncrief’s sound 


metric reduces, up to a constant factor, to the acoustic metric given in Eq. (151. 
For the particular case of the Bondi flow we obtain 


0 = — [—(■Cg — u^)dr — 2udtdr + dr^ + -I- sin^ ddip"^)] . (16) 


The acoustic metric (16) is spherically symmetric and possesses the Killing 
vector field 


^=1 

dt 


(17) 


in the “Newtonian” spacetime M = M x endowed with the metric 0. In 
general relativity, the spacetime metric often provides a natural normalization 
for a timelike Killing vector field k. For example, if the spacetime is asymp¬ 
totically flat, one may require that k has unit norm at infinity. In the present 
case, there is no natural (fluid-independent) spacetime metric to provide the 
normalization; however we may normalize k by requiring k[t] = 1 with t the 
Newtonian time. 

The negative square norm of the Killing vector field k with respect to the 
acoustic metric is given by N{r) := —i3{k,k) = pvs{l — Consequently, 
k is timelike in the subsonic region, spacelike in the supersonic region, and 
null at the sonic sphere r = r^. Since the surface r = Tc is a null surface, 
it constitutes a Killing horizon for the metric (161 ^ The associated “surface 


gravity” k, defined by , will play an important role 

later. It is given by 


K = - V. 


du 

dr 


^oc 

Tr 


5-37 GM 


( 18 ) 


^ See for instance [32II33| for a definition and properties of Killing horizons. 
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In deriving this result, we have differentiated Eq. Q twice and evaluated the 
result at X = Xc in order to compute the first derivative of ly at the critical 
point, which yields 

' - ,19) 


dx 


5-37 


Notice that with our normalization of the Killing vector field and definition 
of the surface gravity, k has units of frequency. This should be contrasted 
to the case of static or stationary black hole spacetimes in general relativity, 
where due to a different normalization involving the speed of light, the surface 
gravity is interpreted as an acceleration, see [32]. For the purpose of the present 
article and the calculation below of the quasi-normal acoustic frequencies, the 


expression (181 not only has the adequate units, but in fact we will see that 


it also provided the correct scaling for the frequencies. 

For the following computations it is convenient to rewrite the acoustic 


metric (16) in diagonal form by means of a suitable coordinate transformation. 


Introducing the new time coordinate 


T = t- 


udr 


the metric outside the sonic horizon assumes the form: 

dr'^ 


e=P- 


-(1 - v^)vidT^ + 


1 — 


+ sin^ ddip^ 


For large radii, a; ^ 7 — 1, we have x'^v ^ \ and p ~ Pao, Vg — 
the acoustic metric becomes (up to a constant conformal factor) 




( 20 ) 


and hence 


-I- -I- sin^ ddip'^) , (21) 


which shows, in particular, that the acoustic metric is asymptotically fiat|^ 
Summarizing this section, we have found that in what the propagation of 


sound waves is concerned, the geometry described by the metric (16) has the 


same qualitative features as a static, spherically symmetric and asymptotically 
flat black hole. The sonic horizon r = r,. plays the role of the event horizon in 
this analogue acoustic black hole. 


4 Quasi-normal acoustic frequencies 

In the previous section, we discussed the geometric properties of the acoustic 
metric for the transonic Bondi flow. Since this metric describes an acoustic 

^ Note that the asymptotic metric j2l[ | is equal to the “canonical acoustic hole” introduced 
in m- Naively, one might think that this metric describes the limit 7 ^ 1 of the metric ( |20[ |; 
however as we will discuss further below the relation = X used in deriving the asymptotic 
form of the metric does not hold for arbitrary values of x in the limit 7 ^ 1 . 















Quasi-normal acoustic oscillations in the transonic Bondi flow 


9 


hole, and since the acoustic perturbations of the transonic Bondi flow are 
governed by a wave equation on the effective geometry described by this metric, 
it is not surprising to find that when perturbed, this flow exhibits quasi-normal 
acoustic oscillations. In the relativistic case, we showed that such oscillations 
arise from fairly arbitrary initial excitations |24j by solving numerically the 


Cauchy problem for the wave equation (14). Here, we do not compute the 


full solution of the wave equation, but rather we focus on the quasi-normal 


mode solutions of Eq. (14) and compute the associated complex frequencies 


s = cr -b iw S C based on our matching method described in 
The quasi-normal modes have the form 

6'ip= 


( 22 ) 


where R := sjplv^r the areal radius associated with the acoustic metric, 
4>{s,r) is a complex-valued function to be determined, and denote the 
standard spherical harmonics with angular momentum number £m. Introduc¬ 


ing the ansatz (221 into Eq. (14) and using the diagonal parametrization (20) 


of the acoustic metric, we obtain the following “eigenvalue” problem: 






-b [s^ +N{r)Vi{r)\ (/) = 0, 


with the function N and the effective potential Vi given by 

= Us(l - v^), Vi{r) = Vs -^ ■ 


R dr 


dr 


(23) 


(24) 


The quasi-normal frequencies are determined by the following requirement on 


4> [23 : Let <j)+{s,r) and 4'-{s,r) be the unique solutions of Eq. (23) satisfying 
the boundary conditions 


lim e®’’*(()+(s, r) = 1, lim e (s, r) = 1, 

r*—>-oo r»—¥ — oo 


r 



for each Re(s) > 0, and consider their analytic continuation to the left complex 
plane Re(s) < 0. Then, the quasi-normal frequencies are determined by those 
values of s for which the functions (j)+{s, ■) and (j)- (s, •) are linearly dependent. 
In [24j we discussed a new numerical method for computing the analytic con¬ 
tinuations of the functions (()±(s, •) based on an iteration scheme, where each 
iteration involves computing numerically an integral over a suitable path in 
the complex r-plane. Based on this method, the quasi-normal frequencies can 
be found by determining the zeros of the Wronskian 


W{s) := det 


(t>+{s,r) (j)-{s,r)\ 


based on a standard Newton algorithm [36] . 














10 


Eliana Chaverra, Olivier Sarbach 


In terms of dimensionless quantities, Eq. (231 can be rewritten as 






+ [s^ + Af{x)Vi{x)] (j) = 0, 


(26) 


with 


Af{x) 

Ve{x) 


{x'^vY 

1 


1 


(x^i^Y 

7-3 


£{£ + 1 ) — 2 ( t (1 — ly^) — a(l + —- 

V 


4(7 + 1) 


(3- 


fYY 

V v 


7-3 ^^2-^ x‘^Vxx 


2(7 + 1) 


(27) 


.(28) 


Here an in the following, and denote the first and second derivatives of 
V with respect to x. Since asymptotically z/ ~ A/a:^, the effective potential Vi 
decays as fast 1/a;^ when a; —)■ oo. In fact, it follows that x^Vi{x) = 

£(£+l)/A^'^, which shows that the centrifugal term dominates for large x. The 
first derivative can be computed from 


XVx 

V 


= A 


- 2 <t xgx{x) 


as long as X Y whereas for x = x^ one can use the expression (19). The 
second derivative can be computed using the identity 

(A X g^Yx) - fyu{v)x Vx) 1 

which follows from differentiating Eq. ([^ twice, and this also shows that the 
effective potential is regular at the sonic horizon. In Eq. (26) the dimensionless 
frequency s is related to s via the relation 


= A” 


GM 


which shows that s must scale like the surface gravity k of the acoustic black 
hole, see Eq. (18). 

For the computation of the functions J\f{x) and Vi{x), one needs to de¬ 
termine as a function of x, which has to be determined by inversion of 
the function / in Eq. ([^. We invert / numerically, using a standard Newton 
algorithm (SS], starting from the seed value 


uq{x) = 


[1 + B+{x - Xc)][l + B^{x - Xx)] 


B+ := 


1± a/1- 


5-37 

2A 


5-37 

2 


which has the correct asymptotic behavior for small \x — Xc\ (see Eq. @) and 
for large a:, for which v ~ Xjx^. 
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Results for the frequencies of the fundamental quasi-normal acoustic os¬ 
cillations of a polytropic fluid with adiabatic index 7 = 1.3333 ~ 4/3 and 
different values for the angular momentum number £ are shown in Table 
Note that for high i the real part of s seems to converge to a value close to 
—0.39k while the imaginary part augments by a value close to 0.6k as £ in¬ 
creases by one. As we will discuss in the next section, these features can be 
understood from simple analytic considerations in the high-frequency limit. 



s/k 

zi(Im(s)/K) 

1 

-0.3863 + 0.8217* 


2 

-0.3887 + 1.455* 

0.633 

3 

-0.3893 + 2.069* 

0.614 

4 

-0.3896 + 2.678* 

0.609 

5 

-0.3898 + 3.284* 

0.606 

6 

-0.3898 + 3.888* 

0.602 

7 

-0.3899 + 4.491* 

0.603 

8 

-0.3899 + 5.094* 

0.603 

9 

-0.3899 + 5.697* 

0.603 

10 

-0.3899 + 6.299* 

0.602 


Table 1 Fundamental quasi-normal frequencies for acoustic perturbations of a Bondi flow 
with adiabatic index 7 = 1.3333 and angular momentum number i. Four significant figures 
are shown. The last column shows the difference between the imaginary part of s/k for the 
given value of I and the same quantity for the previous value of i. 


The frequencies for the first few overtones are given in Table and the 
corresponding spectrum is shown in Fig. 


rio 

II 

s/K {1 = 2) 

s/k {1 = 3) 

II 

s/k {i = 5) 

1 

-1.259 + 0.6782* 

-1.207+ 1.350* 

-1.189 + 1.991* 

-1.182 + 2.615* 

-1.178 + 3.231* 

2 

-2.249 + 0.5629* 

-2.115 + 1.204* 

-2.045 + 1.859* 

-2.001 + 2.503* 

-1.991 + 3.136* 

3 

-3.262 + 0.4965* 

-3.091 + 1.082* 

-2.967+ 1.714* 

-2.888 + 2.362* 

-2.841 + 3.007* 

4 

-4.275 + 0.4539* 

-4.097 + 0.9924* 

-3.938 + 1.587* 

-3.819 + 2.219* 

-3.735 + 2.865* 

5 

-5.287 + 0.4237* 

-5.108 + 0.9279* 

-4.935 + 1.486* 

-4.787 + 2.089* 

-4.670 + 2.722* 

6 

-6.295 + 0.4007* 

-6.122 + 0.8789* 

-5.943 + 1.406* 

-5.777 + 1.981* 

-5.634 + 2.591* 

7 

-6.993 + 0.3942* 

-7.135 + 0.8422* 

-6.955 + 1.342* 

-6.778 + 1.889* 

-6.619 + 2.477* 

8 


-8.144 + 0.8087* 

-7.968 + 1.289* 

-7.787 + 1.814* 

-7.618 + 2.378* 

9 


-9.155 + 0.7825* 

-8.980 + 1.246* 

-8.798 + 1.749* 

-8.622 + 2.292* 


Table 2 First few overtones of the acoustic frequencies shown in the previous table. Here, 
no stands for the overtone number, rio = 1 denoting the first overtone, rio = 2 the second 
one, etc. 


5 Eikonal limit 

As we show in this section, the quasi-normal frequencies found numerically 
in the previous section can be approximated to a fairly good accuracy using 
the high-frequency limit, in which case the quasi-normal frequencies can be 
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g/k 


Fig. 1 The spectrum of the quasi-normal acoustic frequencies of a polytropic Bondi flow 
with adiabatic index 7 = 1.3333. 


obtained directly from the properties of unstable circular null geodesics. As 
discussed in [3I1I3S] and references therein, in the high-frequency limit the 
quasi-normal oscillations can be interpreted in terms of decaying wave packets 
which are localized along an unstable circular null geodesic. Denoting by f^circ 
and Acirc the angular velocity and the Lyapunov exponent, respectively, of the 
unstable null geodesic, the quasi-normal frequencies for high values of € are 
given by the formula m 


s 



Ar 


i£ fir 


(29) 


with rio the overtone number. For an arbitrary asymptotically flat, static spher¬ 
ically symmetric metric of the form 

ds^ = -F{r)dT^ + -t idd^ f sin^ (30) 

G(r) ' ' 

with time coordinate T and positive smooth functions F(r) and G{r), there 
exists an unstable circular null geodesic at r = Vcirc if and only if the function 
F[(r) •= F{r)/r^ has a local maximum at r = Tdrc, and in this case the 
associated angular velocity and Lyapunov exponent are given by m 


Gcirc \/Gd (t circ^ ; Ar^rc ^ 


F{r)G{r) 


1 d2 

H{r) dr'^ 


H{r) 


(31) 


Comparing Eq. (301 with Eq. (20) and discarding the conformal factor p/vs 
which does not affect the null geodesics as trajectories in spacetime we find 
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H{r) = v1{l — yjF{r)G{r) = such that 

^circ = H{Xcirc), 


A - — V 

■^^r.rrr. — 


GM 
vL l-v"^ 


1 d2 


GM [x'^vY y 2H{x) dx'^ 


H{x) 


(32) 

(33) 


with the dimensionless function 

^ia+2^2a ’ 

where v is implicitly determined by Eq. Q. A plot of H{x) for the typical 
value 7 = 4/3 is shown in Fig. and as is visible from this plot there is a 
unique unstable circular null geodesic corresponding to the maximum of this 
function. Computing this maximum numerically, we find, for 7 = 4/3, 



Fig. 2 A plot of the function H{x) for 7 = 4/3. There is a maximum at x = Xdrc — 0.35396 
where u = Udrc — 0.76842. 


Ar 


2k 


~ 0.39006, 


a 


~ 0.60095, 


(34) 


which agrees very well with the corresponding values in Table In fact, the 
formula (29) also fits reasonably well with the first few overtones in Table 
For example, for rio = 3 Eq. (29) predicts Re(s)/H: ~ —2.73042 which agrees 
with the corresponding value for £ = 5 in Table to about 4% accuracy. For 
higher values of rio, a good fit seems to require higher values of 1. 
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6 Dependency on 7 

Finally, we analyze the dependency of the fundamental frequencies on the 
adiabatic index 7 . Table shows such frequencies for £ = 1,2,3 and four 
different values of 7 , varying in the interval 1 < 7 < 5/3. The results from 
the eikonal approximation described in the previous section are also given for 
comparison. As for the case 7 = 4/3, the Lyapunov exponent provides a very 
accurate description for the decay rate, and the angular velocity gives a very 
good description for the increment in the oscillation frequency as I increases 
by one. 


7 

s/k = 1) 

s/re (£ = 2) 

s/re {£ = 3) 

^circl (2/^) 

^circ! ^ 

1.66 

-0.4345 -1- 2.260i 

-0.4343 + 3.924* 

-0.4343 + 5.554* 

0.4343 

1.604 

1.50 

-0.4048 + 0.9858* 

-0.4053 + 1.732* 

-0.4055 + 2.459* 

0.4057 

0.7125 

1.15 

-0.3698 + 0.7335* 

-0.3743 + 1.307* 

-0.3755 + 1.861* 

0.3768 

0.5414 

1.01 

-0.3591 + 0.6905* 

-0.3650 + 1.234* 

-0.3666 + 1.759* 

0.3681 

0.5121 

1.00 

-0.3584 + 0.6879* 

-0.3644 + 1.229* 

-0.3660 + 1.753* 




Table 3 Fundamental quasi-normal acoustic frequencies for different values of 7. Four sig¬ 
nificant figures are shown. The values for 7 = 1.00 are obtained using the limiting procedure 
described below. 


In Fig. 1^ we provide plots for the dependency of the decay rate and oscil¬ 
lation frequencies on 7 . Empirically, we found that 


bJ 

K 




£=1,2,3, 


with Ci a constant depending on I and /3 a constant which is approximately 
equal to 0.25. In particular, this implies that w/k diverges as 7 —)■ 5/3. This 
makes sense since (at the Newtonian level considered here) there is no tran¬ 
sonic, steady-state radial accretion flow when 7 > 5/3. 



y 


Fig. 3 Fundamental quasi-normal acoustic frequency as a function of the adiabatic index 
7 in the range 1 < 7 < 5/3. Left panel: real part of s. Right panel: imaginary part of s. 
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In contrast to this, a /k and w/k have hnite limits when 7 —>■ 1. This limit 
can be understood by first rewriting Eq. 0 in the form 


~Y' 


( _ \ 1 /2cr 


and taking the limit 7 —>■ 1 on both sides, keeping x fixed. Noting that A —)■ 
e^/^/4, and assuming that v converges to a finite value in this limit, we obtain 
the implicit relation 

4 


between x and v. As a consequence, (x^vY 
in Eqs. ( 27|28 ) simplify to 


1 , and the functions N and Vt 


M{x) = I — 


VYx) 


EY + I) 

x"^ 




1 

H-1 

X 


Computing the fundamental quasi-normal frequencies directly with these ex¬ 
pressions yields the results shown in the last line of Table which provide the 
correct limiting values as 7 —> 1 . 


7 Conclusions 

We have analyzed linear acoustic perturbations of the Michel flow solution, 
describing a transonic radial accretion flow into a nonrotating black hole. 
Whereas previous work by the authors [24] treated this scenario using the 
full relativistic Euler equations, in the present article we focused on the re¬ 
alistic case where the sound speed at infinity is much smaller than the speed 
of light. In this limit, the accretion flow at and outside the sonic horizon is 
accurately described by the nonrelativistic Euler equations, and thus a New¬ 
tonian description is sufficient for an understanding of the behavior of the flow 
outside the sonic horizon. The solutions for the unperturbed flow in this case 
have been given by Bondi El, and the dynamics of its acoustic perturbations 
are described in an elegant and efficient way by a wave equation on an effective 
curved background geometry described by the acoustic metric introduced by 
Unruh [28]. For the transonic Bondi flow, the causal structure described by 
the acoustic metric has the same properties as a static, spherically symmetric 
black hole whose event horizon is described by the sonic horizon. Therefore, 
as far as the propagation of sound waves is concerned, the acoustic geometry 
has the same qualitative properties as a black hole, and thus it describes an 
analogue black hole [39|. 

In analogy with a scalar field propagating on a stationary black hole back¬ 
ground, the acoustic perturbations of the Bondi flow exhibit quasi-normal 
oscillations. Using the analogue black hole description, we have computed the 
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associated quasi-normal frequencies s = a + iu, where cr < 0 describes the 
decay rate and lo the oscillation frequency. In the fully relativistic case, we 
had computed these frequencies numerically |24j and found empirically that 
they scaled like the surface gravity k of the acoustic black hole in the limit 
where the sonic horizon is much larger than the Schwarzschild radius of the 
black hole. In the present article we have managed to provide a proof for this 
behavior by reconsidering the problem in the Newtonian limit and by showing 
explicitly that s scales like k. In particular, our Newtonian calculations show 
that for high values of the angular momentum number the fundamental 
quasi-normal acoustic frequencies of a polytropic fluid with adiabatic index 
7 = 1.3333 ~ 4/3, the real part of s converges to a value close to —0.39k 
while the imaginary part augments by a value close to 0.6k as I increases by 
one. As we further showed in this article, these features can be understood 
from an analysis in the high-frequency limit. Therefore, our results provide an 
explanation for the empirical formula (45) given in Ref. [21] . 

We have also analyzed the dependency of the quasi-normal acoustic fre¬ 
quencies on the adiabatic index 7 of the flow. We found that the decay rate and 
oscillation frequency divided by k, ct/k and w/k, increase in magnitude as 7 
increases from 1 to 5/3 and have finite limits as 7 —?► 1. In contrast to this, the 
numerical data suggests that w/k diverges polynomially as 7 approaches the 
value of 5/3, which makes sense since the transonic Newtonian flow considered 
in this article ceases to exist for 7 = 5/3. 

We end this article with brief comments on the feasibility of observing or 
measuring the quasi-normal acoustic oscillations discussed here. First, let us 
comment on the frequency scales associated with these oscillations. As we have 
shown, both the decay rate and the oscillation frequency scale like the surface 
gravity k of the acoustic black hole. As follows from Ref. Plj, the surface 
gravity for a fully relativistic, polytropic flow accreted by a Schwarzschild 
black hole is 


K = 


9 + ( 5 - 37 )( 4 ec- 3 ) 1 c 
32 


(35) 


with = Tc/th the ratio between the sonic horizon radius and the Schwarzschild 
radius th = 2GM/c^. When the sound speed at infinity Voo is much smaller 
than the speed of light c, we have, for 1 < 7 < 5/3, — (5 — 37 )/( 8 u^/c^), 

and the above expression for k reduces to the Newtonian expression given in 
Eq. (18). For 7 = 4/3 in particular, k is of the order of 



with Mq the solar mass and Voc the sound speed at infinity. Comparing 
this expression with the corresponding scale for the fundamental quadrupolar 
quasi-normal frequency of a Schwarzschild black hole (see [38] and references 
therein). 



Ibh ~ 10 


kHz, 
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we see that the fluid oscillations’s frequencies are suppressed by a factor of 
about 100(uoo/c)^ compared to the quasi-normal frequencies associated with 
the black hole. For the situation where gas is accreted from the interstellar 
medium, where Vao — lO/cm/s, this factor is extremely small, leading to very 
large periods for the fluid oscillations even for stellar-mass black holes. The 
situation improves somehow when considering a fluid flow with adiabatic index 
7 = 5/3, in which case — 3/(8uoo/c) when Voo ^ c, see for instance [23]. In 
this case, we obtain 


8^2 c 

3 th 




kHz, 


which is quadratic instead of cubic in v^ajc', however, for gas accreted from 
the interstellar medium this is still very small compared to the fundamental 
quadrupolar quasi-normal frequency of the accreting black hole. 

Second, when the backreaction of the fluid on the metric field is taken 
into account, the quasi-normal acoustic oscillations with angular momentum 
number £ > 2 should produce gravitational radiation which could be mod¬ 
eled, to a good approximation, by considering linear metric perturbations of 
a Schwarzschild black hole (see Ref. [10] and references therein for an efficient 
formalism for describing such perturbations). For a related recent study on 
gravitational waves triggered by a non-spherical dust cloud which is accreted 
by a Schwarzschild black hole, see mi- However, in the scenario discussed in 
this article it does not seem that the amount of gravitational radiation will 
be significant since (i) we only consider small departures of a spherically sym¬ 
metric flow and (ii) because the characteristic frequencies associated with the 
acoustic perturbations are so much smaller than the characteristic frequencies 
of the accreting black hole. On the other hand, the situation could change 
dramatically when considering more realistic flows like those occurring in ac¬ 
cretion disks, where the flow has already nonvanishing angular momentum in 
the unperturbed solution. We intend to generalize our work to such scenarios 
in future work. 
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